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Evaluation of Exponentially Weighted Semi-Infinite Integrals 
by Quadrature (Laguerre Quadrature) 
Closed subroutine , with one program parameter 
18 + 2N (see below) 

N(1.8 + T). milliseconds, where T is the toation in milli- 
seconds of the auxiliary subroutine. 
Location (may be used by auxiliary subroutine) 
When this routine is located at y , entry is made by the orders 
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where a is the location of the auxiliary subroutine which 
computes the values of the function to be integrated. When 
control is returned to the right side of p + 1, the computed 
integral will be in the accumulator register and location 

y + 1^. 

To evaluate the integral 

CD 

e' X f (x) dx, 
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this routine uses a form of Gaussian Quadrature appropriate 
to the interval (0,00) and the weighting function e 

J" e- X f(x)dx~ ^^f^). 



(1) 



The values A. and i are chosen in a manner such as to give 
no truncation error when f (x) is a polynomial of degree 2N -1 



or less. In the case where the factor e 
explicitly in the integrand, 



-x 



does not occur 
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CD -oo t N 



g(x) dx = J e" x [e x g(r)] dx Z L. Z ^ e K g^) 



It is assumed that the function e g(x) may be closely- 
approximated by a polynomial function. 

Because the actual values of the points X, and the weights 
A, and B may exceed 1, they have been scaled down by powers 
of two# P and Q are defined in equations (l) and (2), and 
R is defined below. 
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The auxiliary subroutine which computes f (x, I must take 
the scaling of these values of x, into account. The 
function values computed by the auxiliary are assumed to 
lie in the range -1 < f (x, ) < 1. 

The closed auxiliary subroutine is entered from the main 
routine with x, in the accumulator and link in Q; control 
is returned to the main routine with f (x, ) in the 
accumulator. 
USE: To use this routine, the programmer copies the Integration 

routine first on his program tape, and immed lately after 



SCALING: 



- 3 - 

it the parameters, points x fc , and weights ^ or B k appropri- 
ate to his needs. These latter numbers appear on the tail 
of the library tape, labeled by the number N of points at 
which the function is to be evaluated, and the type of weights 
(A, or B, ) to be used. 

The scaling of the values of x, is such that the auxiliary 
subroutine is presented with x, *, where 

x k * = 2" R x k , 1 < R < 6, 

and the largest x, satisfies 1/2 < (x, ) max < 1 for all N. 

P/ Q\ 

The computed integral is scaled down by 2 (or 2 ;. 



The A, for N = 9 to 15 are not included because they decrease 
rapidly with increasing k and N, and become too small to 
be held in a single Illiac register. 

For the convenience of the programmer, the above scale factors 
are contained in the subroutine parameter at location y+l6, 
in the following form: 



(y+16) = OR (y+18) OP (y+l8+N), 
= OR (y+18) OQ (y+18+N) 



for A weights 
for B weights. 



ACCURACY: 



REFERENCES : 



The truncation error due to omission of powers of x higher 
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There will also be round-off errors which may be significant; 
for a discussion of these see the write-up of library routine 
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Preset point and weight addresses 
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Plant subroutine entry, and link 
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clear sum box 
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Save L.S.P. of sum 
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get x,; link to Q 
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jump to auxiliary 
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L.S.P. to A, weight in Q 
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x f(x k ), accumulate M*S. P. 
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step x, address 
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step A, address 
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exit via link 
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temporary store for M.S. P. 
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Ceofficient addresses and scale constants 
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end constant 



